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Properties of Lamb Waves Relevant to the Ultrasonic 
Inspection of Thin Plates* 


T. N. GRIGSBY t, MEMBER, IRE, AND E. J. TAJCHMANT 


Summary—An experiment is described, using 1/32-inch thick 
steel plates, in which the insertion loss of artificial flaws was 
measured. The artificial flaws consisted of 1/32-inch wide saw cuts 
of depths varying from 0.000 inch to 0.021 inch in increments of 
0.003 inch. No simple relationship was found relating the insertion 
loss to the depth of the saw cuts, although the presence of a flaw 
was usually indicated. Measurements of the group velocity of Lamb 
waves are described and the values obtained are shown to agree 
with calculated values. The theory of Lamb waves is discussed, 
and the results of computation of the frequency equation for the 
symmetrical case for a ratio of longitudinal velocity to shear ve- 
locity of 1.8 are given in the form of curves of phase velocity vs 
frequency-thickness product. Curves of group velocity vs frequency- 
thickness product for this case are also given. A formula is derived 
which gives group velocity as a function of phase velocity and 
frequency-thickness product. A derivation of the frequency equation 
is given. 


INTRODUCTION 


OR the purpose of nondestructive testing of thin 
Frites with thickness of the same order of magnitude 

as the wavelength, it is desirable that as much 
knowledge as possible be acquired about the possible 
modes of ultrasonic propagation in such plates. Although 
techniques of using Lamb waves for ultrasonic inspection 
were patented in 1951,' and although they have been 
applied in some cases,’ they have failed to gain wide- 
spread acceptance. It is the purpose of this paper to 
describe some of the properties of Lamb wave propa- 
gation which are relevant to possible nondestructive 
testing applications. 

The case of waves of ultrasound being guided by a 
thin elastic plate is, in general, qualitatively analogous 
to the case of electromagnetic waves being guided by a 
rectangular wave guide. In the case of ultrasound being 
guided by a thin plate, the modes of propagation are 
solutions to the equations of motion subject to the bound- 
ary conditions imposed by the medium in contact with 
the faces of the plate. In the corresponding case of electro- 
magnetic waves being guided by a rectangular wave- 
guide, the modes of propagation are solutions of Maxwell’s 
equations subject to the boundary conditions imposed 
by the walls of the guide. Lamb waves include only those 
modes of vibration of an infinite plate in which the 
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particle motion is independent of the direction which is 
parallel to the free surfaces and normal to the direction 
of propagation. However, in addition to Lamb waves, 
there are other modes which have particle motion only in 
the direction which is parallel to the free surfaces and 
normal to the direction of propagation. For the first and 
higher modes of this latter type, the frequency equation 
is identical to the frequency equation for the propagation 
of electromagnetic waves in a rectangular wave guide.’ 
The nondispersive zeroth order mode of this case docs not 
have a counterpart in the rectangular waveguide.* No 
attempt has been made so far to use these modes in non- 
destructive testing. 

From the standpoint of nondestructive testing, the 
main features that distinguish Lamb waves from longi- 
tudinal and shear waves are that Lamb waves are dis- 
persive and that there is an infinite number of possible 
modes. 


Lams WAVE PROPAGATION 


Lamb waves are waves which can exist in an isotropic, 
lossless, elastic plate of infinite length and width, but of 
finite thickness, bounded on both of its faces by vacuum 
and which are independent of a direction parallel to the 
free surfaces and normal to the direction of propagation. 

Several theoretical treatments of Lamb waves may be 
found in the literature." The frequency equation relat- 
ing phase velocity to frequency can be obtained by 
making use of Hooke’s Law, Newton’s second law of 
motion, and by satisfying the boundary conditions that 
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the tensional stress normal to the free surfaces and the 
shear stress tangential to the free surfaces must vanish. 
It can be shown that the frequency equation for Lamb 
waves may be split into two separate equations one of 
which applies to motions of the plate which are sym- 
metrical about the medial plane and the other to motions 
which are anti-symmetrical about the medial plane.” 
For the convenience of the reader, this derivation appears 
in Appendix I. 


The frequency equation for the symmetrical case is’? 








san (2) 
ies E _ e+ By. a 
(2) th a8 i 
tanh D, 
where 
T? 
a = kafl — Te? 
ER E 
in which 
f = frequency 
h = thickness of plate 


V = Lamb wave phase velocity 


k = wave number = 2rf/V 

V: = velocity of longitudinal waves in an extended 
medium of the same elasticity and density as 
the plate 

Vs = velocity of shear waves in an extended medium 


of the same elasticity and density as the plate. 


The frequency equation for the antisymmetrical case is!* 


tanh (24) A 
ees Ne a 2 
beak (2 ) (k? +B")? 


2 

The frequency equations are solved by assuming a 
value of V and then computing the corresponding values 
of frequency-thickness product by successive approxi- 
mations, preferably with the use of a digital computor. 
Worlton has found roots of these frequency equations for 
zirconium, stainless steel, uranium, brass and aluminum. 
and has plotted the results as curves of phase velocity vs 
the frequency-thickness product.’ Firestone and Ling 
previously made calculations for aluminum and also gave 
curves of the group velocity. Normalization of the fre- 
quency equation permits the computation of a single set 
of curves for each of several values of the ratio of longi- 
tudinal velocity to shear velocity so that the desired 
solution to the frequency equation could readily be found 
from such curves for plates of any isotropic solid. A set 


2 H., Lamb, op. cit., p. 116. 

133 Tbid., p. 122. 

4D. C. Worlton, “Lamb Waves at Ultrasonic Frequencies,” 
Hanford Atomic Products Operation, Richland, Wash., Rept. No. 
HW-60662, pp. 31-32, pp. 51-58; 1959. 
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of normalized curves showing phase velocity vs frequency- 
thickness product for a ratio of longitudinal velocity to 
shear velocity of 1.8 (i.e., Poisson ratio = .277) is shown 
in Fig. 1. 

It is easily shown (Appendix IT) that the group velocity 
may be expressed in terms of phase velocity and frequency- 
thickness product, as follows: 

dV 


‘ 
R= T 
d{fh) 


To normalize this expression with respect to shear wave 
velocity, let Von = V,/Vs, Vn = V/V s and (fh), = fh/Vs; 
then, 


1 
1 — (fh) 
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Vin = Fay baee |. + 
i area 


which is the formula used to compute the curves of Fig. 2. 
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Fig. 1—Normalized curves of phase velocity vs frequency-thick- 
ness product for a ratio of longitudinal velocity to shear velocity 
of 1.8. 
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Fig. 2—Normalized curves of group velocity vs frequency-thick - 
ness product for a ratio of longitudinal velocity to shear velocity 
of 1.8. 


EXPERIMENTAL APPARATUS AND MATERIALS 


A functional block diagram of the apparatus is shown 
in Fig. 3. The following equipment was used: gate pulse 
generator, Electropulse 4120A; pulsed oscillator, Arenberg 
PG-650-C; transducers, Branson type Z, 0.5 Me, 1.6 Me, 
3.3 Mc, and 10 Me, 0.5 inch in diameter; preamplifier, two 
Hewlett-Packard type 460A wide-band amplifiers in 
cascade; attenuator, Arenberg ATT 693; amplifier, 
Tektronix No. 121; oscilloscope, Tektronix 545 with 
53/54C dual-trace preamplifier. The wavemeter was 
constructed using quartz crystals resonant at 0.5000 Me, 
1.6000 Me, 3.3000 Me, and 10.000 Me. 

The test plates consisted of precision-ground flat stock, 
Starrett type No. 496. Fig. 4 shows a 10 X 18-inch test 
plate; however, most of the data given below were obtained 
using 2 X 10-inch test plates, since no side reflections were 
observed for this size of plate. The test plates were 
supported in an aluminum frame so that they were 
bounded on both faces by air. 

The positioning apparatus and lucite wedges used are 
shown in Fig. 4. 

In order to generate Lamb waves of phase velocity, V, 
one can project a longitudinal wave onto the plate at an 
angle of incidence of 


6 = sin™ (V/V), 


where V, is the velocity of the longitudinal wave." 
Tor the data which follow, FV, velocity of longitudinal 
waves in lucite = 2680 m/sec. 

The couplant used was ordinary motor oil. 


MEASUREMENT OF GROUP VELOCITY 


The method used for measuring group velocity was to 
set the transmitting and receiving transducers close 
together and measure the distanee and delay time. The 
transmitting and receiving transducers were then moved 
further apart, and the distance and delay time were again 
measured. The difference between the two distances 
divided by the difference between the two delay times is 
the group velocity. There is an inherent uncertainty in 
the measurement of group velocity. Since Lamb waves 
are dispersive, the different components of the pulse travel 
with different velocities. Consequently, the transmitted 
and received pulses do not have the same shape (Fig. 4), 
and some arbitrary criterion must be established for 
deciding at what point in time a given pulse occurs. When 
the data of Table I were gathered, the time at which a 
given pulse occurred was taken to be the time at which 
its leading edge reached 50 per cent of the maximum 
voltage attained. The agreement between the calculated 
and experimental values of group velocity is good when 
one considers that the ratio of longitudinal to shear 
velocity in steel is actually 1.85 instead of 1.8, which was 
used in computing Fig. 2. 


15 D. C. Worlton, ‘Ultrasonic testing with Lamb waves,” Non- 
destructive Testing, vol. 15, p. 219; July-August, 1957. 
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Fig. 3—Functional block diagram of apparatus. 








Fig. 5—Input waveform (lower trace, actually the output of the 
wavemceter of Fig. 3) and output waveform (upper trace) resulting 
from excitation of the second asymmetrical mode in a 1/32-ineh 
thick steel plate at 2.46 Me. The pulse duration is 15 ysec, and 
the total sweep time is 200 usec. 


TABLE I 























Plate Group velocity, ips 
thick- Fre- Phase 
ness, quency, velocity, Experi- Calculated 
inches Me. ips mental (V,/Vs5 = 1.8) 
| 1 
1/8 3.44 308 , 500 75,000 70,500 
1/32 2.99 164,100 69,900 71,200 
1/32 9.13 164, 100 90,650 | 88 , 900 
1/32 4.95 228,000 188,750 195,200 
1/32 5.40 228 , 000 190, 750 190,000 
1/32 8.70 233 , 800 180 , 900 19+, 000 
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MEASUREMENT OF INSERTION Loss 
oF IpEALIzED FLAWS 


In order to obtain data which might indicate possible 
applications of ultrasonic Lamb waves in the nondestruc- 
tive testing of thin plates, an experiment was conducted 
using plates having highly idealized flaws. The primary 
objectives of this experiment were to determine whether 
or not such flaws would have significant insertion losses, 
and, if so, to attempt the correlation of insertion loss with 
flaw size. 

The plates used in this experiment were 2 X 10 X .033- 
inches thick. The simulated flaws consisted of various 
depths of saw cuts, 1/32-inch wide. A typical plate with 
its saw cut is shown in Fig. 6. Nine such test plates were 
prepared from one 18 X 10 X 1/32-inch steel plate. Saw 
cuts which varied in depth from 0.003- to .021-inch in in- 
crements of 0.003 inch were made in seven of these plates. 

The experiment performed with these plates consisted 
of comparing the transmission of Lamb waves in plates 
without saw cuts with the transmission in plates with 
saw cuts. The power level was read from the precision 
attenuator when it was set so that the signal at the input 
to the oscilloscope was 0.18 volts peak-to-peak. (If the 
pulse shape was distorted so that the top of the pulse was 
not flat, the 0.18 volts were measured across the average 
envelope of the pulse.) In order to minimize the experi- 
mental error, a standard technique was adopted for 
obtaining acoustic coupling from the lucite wedges to the 
plate. This procedure consisted of applying a spot of oil 
the size of a penny beneath cach wedge, applying enough 
pressure so that some oil was squeezed from beneath each 
wedge and wiping off the excess in front of the wedges. 
The standard deviation of the experimental error was 3 
db when measurements were being made on unslotted 











Fig. 6—Test plate with simulated flaw. 
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reference plates. Consequently, all measurements were 
repeated ten times and averaged so that the standard 
deviation of sample averages is approximately 1 db. In 
some instances, there was considerable distortion of the 
pulse which resulted in some additional error. The results 
of these measurements are presented in Table II. 

The angle given in Table II is the angle of incidence in 
lucite of the longitudinal wave used for excitation. The 
transmitting and receiving wedges were as nearly identical 
as it was possible to make them. 

It is apparent from Table II that no simple relationship 
exists between the depth of the saw cut and the amount 
of insertion loss. It appears that ultrasonic inspection 
with Lamb waves is capable of indicating if there is a 
difference between a given plate and a reference plate. 
It does not appear that Lamb wave inspection will permit 
an unambiguous description of this difference. 

The authors are not able to give a quantitative expla- 
nation of the results of these saw-cut measurements at 
this time. A qualitative explanation is as follows: The 
energy of the wave incident upon the saw cut is con- 
centrated in a single mode which is determined by the 


‘frequency-thickness product of the plate and the trans- 


mitting wedge used. When the incident wave reaches the 
first thickness discontinuity, mode conversion will, in 
general, take place so that all boundary conditions at the 
discontinuity will be satisfied. When that part of the 
wave which is transmitted beneath the saw cut reaches 
the other side, additional mode conversion must take 
place so that all boundary conditions are satisfied. If, 
as was true in these measurements, the receiving wedge 
is similar to the transmitting wedge, the receiver will be 
most sensitive to that mode which was originally trans- 
mitted. The energy at the receiving wedge which is 
effective in producing an output voltage will be approxi- 
mately equal to the energy which was imparted to the 
plate from the transmitting wedge, less the energy which 
was reflected back down the plate in various modes, less 
the energy which arrives at the receiving wedge in modes 
having different phase velocities than the transmitted 
mode. Since the length of the pulse in the plate is much 
greater than the width of the saw cut, a standing wave is 











TABLE II 
Insertion Loss IN DecrBELS oF Various Deprus of Saw Cuts (Saw Cuts Factna Up) 
i 
Frequeney, Me 2.93 9.15 10.7 8.93 9.55 9.50 10.73 
angle 40° 40° 7.5° 20° 30° 20° 15° 
mode 1-S 2-S 3-S 4-S-A 3-A 4-A 5-A 
Depth of l 
cut, inches : 
0.003 0.6 6.8 3.9 12.6 i 5.2 0.6 7.7 
0.006 0.1 15.7 4.9 16.1 | 8.5 1.9 5.7 
0.009 —1.9 2.9 —0.6 10.6 5.4 18.4 * 
0.012 4.5 12.4 4.0 22.9 23.9 16.9 * 
0.015 6.0 4.9 9.9 10.8 10.5 13.1 9.5 
0.018 6.9 13.1 6.8 9.2 5.4 * 8.1 
0.021 3.8 12.8 5.8 20.7 11.5 z 14.5 








*The pulse was too distorted t> measure: 
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present in the plate beneath the saw cut, and its width is 
important in determining its insertion loss. 

The effect of the orientation of the saw cut on its inser- 
tion loss was measured in the following manner: An un- 
slotted reference plate was inserted in the apparatus and 
the frequency was adjusted for maximum amplitude. A 
slotted plate was then placed in the apparatus with its 
slot up and the amount of attenuation required to make 
the signal level at the input of the oscilloscope equal to 
0.18 volts was recorded. The same plate was then turned 
over and the attenuation adjustment was again made and 
recorded. Ten measurements were made and averaged. 
The difference between these two values of attenuation is 
given in Table III. 


TABLE III 


EFFECT OF THE ORIENTATION OF Test PLATES UPON THE TRANS- 
mission Loss or ULTRASOUND—THE TRANSMISSION Loss 
WITH THE Saw Cut Down REFERRED TO THE TRANSMISSION 
Loss WITH THE Saw CUT ur—i1n DECIBELS 





Frequency, Me. 9.4 10.3 9.3 9.15 
angle 40° 27.5° 30° 20° 
mode 2-5 3-5 3-A 4-S-A 

Depth of 

cut, inches 
0.003 0.6 0.2 0.7 6.3 
0.006 1.8 -1.2 0.5 4.0 
0.009 2.0 —0.1 -1.7 — 5.9 
0.012 0.6 —0.9 0.7 —13.8 
0.015 0.0 0.9 —0.7 1.1 
0.018 2.3 1.3 1.6 7.8 
0.021 0.9 —0.4 1.1 — 2.6 














Differences too great to be attributed to experimental 
error occurred only in the case of the composite fourth 
symmetrical and asymmetrical mode. Apparently, more 
of the energy is carried on one side of the plate than on 
the other because of cancellation of the motion on one 
side and reinforcement on the other. 

The frequencies at which the various modes were found 
when the data for Table III were taken differed slightly 
from the corresponding frequencies for the data of Table 
II. This was attributed to small differences in the adjust- 
ment of the angles of the positioning clamps. 


CONCLUSION 


In order for a complete set of curves to be computed 
giving the solutions to the period equation of Lamb 
waves, it is desirable that the period equation be normal- 
ized with respect to shear wave velocity. It would also 
be desirable to have curves of group velocity vs frequency- 
thickness product which are normalized with respect to 
shear wave velocity. 

Measurements were made of the group velocity of 
Lamb waves, and the calculated and experimental values 
agreed satisfactorily. 

Measurements made of the insertion loss of saw cuts in 
0.033-inch thick steel plates indicate that no simple 
relationship exists between the depth of the saw cut and 
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the amount. of insertion loss. This agrees with Intuitive 
concepts of reflection at a thickness discontinuity. Mode 
conversion, which is known to exist at such discontinui- 
ties," is determined by the necessity for satisfying 
boundary conditions at the discontinuity. The analogous 
electromagnetic problem has been discussed by Slater.” 

Measurements were made of the insertion loss of the 
saw cuts when they were on the same side of the plate as 
the transducers and when they were on the side of the 
plate opposite the transducers. A significant difference 
was found in a case where the fourth symmetrical and 
asymmetrical modes were being propagated together. 
This was probably due to cancellation of particle dis- 
placement on one side of the plate and reinforcement of 
particle displacement on the other side so that more power 
was transmitted along one side than along the other. 


APPENDIX I 


DERIVATION OF THE FREQUENCY EQUATION 
FoR LAMB WAVES 


By utilizing Hooke’s law and Newton’s second law of 


‘ motion it is possible to derive the equations of motion for 


a lossless isotropic medium:” 


d'u OA 2 

PB = (A+ n) ag TRV Us 
av ðA é 

pzp =A +u) gy PYT (3) 
aw _ ðA 2 

p aE = Atag + Vw, 


where 
A = Si + Se + S3, 
Ë Ë æ 
= Ox? + ay” + oz. 

In this two-dimensional problem, the x axis is assumed 
to be horizontal with its positive direction to the right, 
and the y axis is assumed to be vertical with its positive 
direction upward. The plate is assumed to be located so 
that it is symmetrical about the zz plane. The type of 
wave motion to be considered is independent of z, że., all 
cross sections of the plate parallel to the xy plane are 
identical, and there is no component of particle motion 
parallel to the z axis. Consequently, it is required to 
satisfy only the first two of equations (3). 


16 F, A. Firestone and D. S. Ling, “Report on the Propagation 
of Waves in Plates—Lamb and Rayleigh Waves,” Sperry Products 
Inc., Danbury, Conn., Tech. Rept. 50-6001, p. 12; 1945. 

1? D, C. Worlton, “Lamb Waves at Ultrasonic Frequencies,” 
Hanford Atomic Products Operation, Richland, Wash., Rept. No. 
HW-G0662, pp. 29-30; 1959. 

18 J. C. Slater, “Microwave Transmission,” McGraw-Hill Book 
Co., Ine., New York, N. Y., pp. 168-173; 1942. 

19 W, P. Mason, “Physical Acoustics and the Properties of 
Solids,” D. Van Nostrand Co., Ine., Princeton, N. J., pp. 12-14; 
1958. 
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It is readily verified that the solution to the equations 
of motion for this two-dimensional case is 


(4) 


where @ and y¥ are auxiliary functions satisfying the 
equations - 
S = (+ WV, 





or 


vy 
P È 


pap =uviv. 6) 

Since 4/(A + 2u)/p = Vi = the velocity of a longi- 
tudinal wave in an extended medium and Vs = V u/p = 
the velocity of a shear wave in an extended medium, (5) 


may be written as follows: 





ay 


— TV2W?2 
ap = VsV¥- 


(6) 


It is assumed that there is periodicity with respect to 
both z and ż; so that ¢ and y may be written, 


$ pe Faje TR, y = GUEST, 7) 
It is readily seen that 
GK F 
T$ = o, (8) 
a 3 
A — k$; (9) 
so that from (6), 
d$ _ (r E sa) 
y? = k Vi $, 
defining 
a= ye on ce ; 
: (10) 
ap > 
y = ag 
Since w = Vk, 
ee 
a= k Ni = re 


Thus, for a given value of phase velocity and wave- 
length, the manner in which ¢ varies with respect to y 
consists of a solution to an ordinary differential equation 
in y. 


Similarly, 
ey = 2 , — . he 
au 8y, where B6=k q|! — FE ay 


The solutions to (10) and (11) are familiar; so that 
Poy = 
Gon = 


ae*” + be*", 


ce” + de®", 


(12) 
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Fn and Gu may be written as hyberbolic functions so 
that 


i(wt—kzr) 


$ = (Asinh ay + Beosh aye 
ý = (Csinh By + Deosh pyjet, 


(13) 


To obtain the frequency equation, it is necessary to 
satisfy the conditions that T, and T, vanish at the two 
values of y = + h/2 where A is the thickness of the plate. 
T, and Ts are given by” 

Ta = MA + QS,, and T, = use. 


A = V’¢ so that from (9) and (10), 








A = (a? — k’*)¢. 
g = &_ ob öy TEn ay 
? ay ay? Ox Oy dx dy 
T, = [\@’ — kA sinh ay + AC’ — k’)B cosh ay 


+ 2uo’A sinh ay + 2ua B cosh ay + 2ikCByu cosh By 
+ 2ikDBp sinh pyle“. 
Collecting terms and noting that 
=Ak? + Aa? + Qua? = uk? + 6°), 
we have l 
Tı = [An(h? + 6”) sinh ay + Bulk? + 8°) cosh ay 
+ C(2ikBy) cosh By + D(2ikBy) sinh Bye?! *?. 
Te = use 


ov, du 8 ee, Ho 
Ox ` dy dxdy ðr” ` dx dy 


[A(—2tkay) cosh ay + B(—2thap) sinh ay 
+ Cu(k? + 6’) sinh By 
+ Dulk? + B) cosh Byje'?'**’. 


(14) 


Ss Ou ay 


ll 


Te 


To make the manipulation less cumbersome, let 
all? + 6°) 


2ikßbu = n, 


sinh (ah/2) 


Il 


—2ikap = p, 
sinh (Bh/2) = sz, 
cosh (Bh/2) 


4, 


Ul 


Say 


cosh (ah/2) = Ca, 


Cg. 


Then, since at y = + A/2, Ta = T = 0, four simul- 
taneous equations in the four unknowns A, B, C, and D 
are obtained?’ 


qa Gla NCp nsa || A 0 

Peas Pa, B PLIO ig 
— Ba (a neg ~ns || C 0 

PCa —PSa — qsg qes ILD 0_ 


20 W., M. Ewing, et. al., op. cit., p. 282. 
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If B and D and the second and fourth columns of the 
square matrix are interchanged, the system may be 
manipulated so that it appears as follows: 


~ 


Ba Ny O 0 A 0 | 
Pea Ca O 0 D 0 


0 0 AC3 Glo 
aa lo 


If this system is to have a nontrivial solution, it is 
necessary that the system determinant vanish, 


(17) 


O O ps Psa 


| Wa NS~p |_| NCg Ca 


= 0. (18) 








Pla ate | Sg PSa 


The system determinant will therefore be equal to zero 
if either of these second-order determinants is zero; 
furthermore, it is clear that any solution of either of the 
second-order systems, 


NS, | | _ | is ar _ 10] 
f , ’ 
Pea w D 0 qe i B lo 
is a solution of the fourth-order system and furthermore 
that all combinations of separate solutions to these two 
second order systems are solutions of the original fourth 
order system. Therefore, it is possible to separate the 
original fourth-order system into two independent second- 
order systems. The two systems are analyzed by assuming 
first that A = = 0 and deriving the resulting period 
equations for the remaining second-order system. Thus, 


$= B(cosh aye’ et? 





(19) 











(20) 
y = C(sinh Bye“, 


which implies motion which is symmetrical with respect 
to the medial plane of the plate, and consequently these 
vibrations are referred to as the symmetrical modes; 
whereas, if it is assumed that B = C = 0 and the other 
period equation is solved the functions are 


= A(sinh aye’, 
y = D(cosh Bye“"'*™, 


(21) 


Since these functions define motions which are anti- 
symmetrical about the medial plane, such motions are 
referred to as the anti-symmetrical modes. 

The frequency equation for the symmetrical case is 
accordingly 


Np Gla | _ 0 
zi $ 








qs PSa 
or 


2 
NPC Sa = Q Cah, 


CaSa 


2 
an A 
Casg np 


March 


Substituting the original expressions baek into the fre- 
quency equation, we obtain 


tanh (ah/2) _ (k° 
tanh (Bh: 2) 


+ BY 
shag 





Similarly, the frequency equation for the anti-sym- 
metrical case is 


tanh (ah/2) 
tanh (Bh/2) 


Akab E 
(k? + By 





(23) 


APPENDIX IT 


DERIVATION oF A FORMULA FoR Group VELOCITY 
IN TERMS oF PHASE VELOCITY 
AND FREQUENCY THICKNESS PRODUCT 


It may be shown that the group velocity is equal to the 
derivative of angular frequency with respect to wave 
number.” 


m Eae 2. 
V, dk (24) 


V + kdV/dk, and since dV /df = 


Since w = Vk, V, 
(dV /dk) - (dk/df), 











ihe 

adv df B -V 

dk kldV v| K] 

E ic "| 

df 
Therefore, 
K 1 | 
Van es eG (25) 

Boa 
a CET 


To normalize this expression with respect to shear wave 
velocity, let V,. = V,/Vs, and since V, = V/Vs and 





(fh), = fh/Vs, (25) may be written 
E eee enn, en 
Vin = Val V, tN 
1— dV, | 
L Uh TD 





This is the equation used to compute the normalized 
curves of group velocity vs frequency-thickness product 
shown in Fig. 2. 


2 S, Goldman, “Frequency Analysis, Modulation and Noise,” 
McGraw-Hill Book Co., Inc, New York, N. Y., p. 117; 1948. 
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APPENDIX III 


List or NOTATION 


1, 0, ¢, d = constants 


A, B,C, D = constants 


a 


Un 
B 


Ss 


cosh (ah/2) 
cosh (8h/2) 
frequency 


normalized frequency-thickness product 


amplitude of @ 
amplitude of y 
plate thickness 
complex operator 
V1 

wave number 


2rf 
V 


2ikßu 





— 2ikau 
a(k? + 6°) 
tensional strain along the x axis 


ðu 
ax 


tensional strain along the y axis 


av 
ay 


tensional strain along the z axis 


ðw 
Oz 


shearing strain about the z axis 


ðv 
Ox 


sinh (ah/2) 
sinh (6h/2) 


du 
dy 


time 
tensional stress parallel to the y axis 


shear stress about the z axis and parallel to the x 
axis 


u 


>= 
li 


<e © d E 
ll 


€ 
Ii 


33 


displacement parallel to the z axis of an elemental 
cube of volume dx dy dz 


displacement parallel to the y axis of an elemental 
cube of volume dz dy dz 


Lamb wave phase velocity 
velocity of an incident longitudinal wave 


Lamb wave group velocity 


= normalized group velocity 


V,/ Vs 


velocity of longitudinal waves in an extended 
medium of the same elasticity and density as 
the plate. 


V/Vs = normalized phase velocity 


velocity of shear waves in an extended medium 
of the same elasticity and density as the plate 


displacement parallel to the z axis of an elemental 
cube of volume dz dy dz 


coordinate in the direction of propagation 
coordinate normal to the faces of the plate 


coordinate parallel to the faces of the plate and 
normal to the direction of propagation 


dilation = S, + So + S; 


angle of incidenee, measured to the normal to 
the faces of the plate, of an incident longitudinal 
wave 


Lamé’s constant 

shear modulus 

density 

auxiliary function for dilation 
auxiliary function for rotation 


angular frequency 


ð“ ð` om 
ant oy a 
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